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The problem of direct interband transitions in semiconducting crystals 
is considered in detail. The solution of the problem is shown to lead to a 
system of simultaneous coupled nonlinear differential equations. These 
equations do not distinguish between static and dynamic fields, and are 
obtained without recourse to any approximation procedures. They therefore 
apply for electromagnetic fields of arbitrary amplitude. It is shown, finally, 
that the linearization of these equations correctly predicts the fundamental 
optical absorption edge and the photoconductive rate equations which result 
in photomixing phenomena. 

I. INTRODUCTION 

Within the past year there has been an increasing interest in semi- 
conductor behavior, resulting from the development of the optical maser. 
In addition to the fact that the semiconductor can be used as an optically 
active medium, this interest results, in part, from a variety of nonlinear 
effects which can occur in semiconductors at infrared and optical fre- 
quencies. 

Optical masers are capable of providing very high-intensity fields at 
photon energies corresponding to the energy gap between valence and 
conduction bands. Since nonlinear effects generally vary as some power 
of field strength, the laser has stimulated keen interest in semiconductor 
nonlinear phenomena. Examples of such phenomena are the Franz- 
Keldysh 12 effect and the multiple-photon process. 3 Closely allied with 
these nonlinear interband effects is the process of photomixing 45 (photo- 
conductivity), where two coherent optical signals are beat together to 
produce a photosignal which is proportional to the instantaneous density 
of photoexcited electrons and hence contains sum and difference fre- 
quency terms. 

The aforementioned phenomena have been analyzed previously by 
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calculating the interband transition probabilities from time-dependent 
perturbation theory 1 ' 3 for the case where an electrostatic and electro- 
magnetic field act on an electron in a crystal. This approach involves 
discarding quadratic terms in the vector potential in the Hamiltonian 
(the static and dynamic fields are treated differently), and neglecting 
higher-order nonlinearities in the population distribution. 

In this article, a more general theory of direct transitions in semi- 
conductors is presented. A semiclassical approach is used, wherein the 
semiconductor properties are quantized and the electromagnetic field 
is treated classically. A set of rather simple coupled nonlinear differential 
equations is obtained for the equations of motion governing the funda- 
mental absorption process, without resorting to any perturbation or 
approximation procedures. In contrast to the usual method of analysis, 
these equations do not distinguish between static and dynamic fields. 
The linearization of the nonlinear equations correctly predicts the funda- 
mental optical absorption, and is shown to lead to the photoconductive 
rate equation for the generation of conduction band electrons. 

II. SOLUTION TO THE WAVE EQUATION 

The equation of motion for an electron in a crystal is the wave equation 

3C^ = ih(df/dt), (1) 

where 3C is the Hamiltonian operator and \f/ is the one-electron wave 
function. Without the presence of the radiation or static fields the 
Hamiltonian, 3Co , is 

3Co = (p 2 /2m) + V(r). 

Here m is the electron mass, V(r) is the lattice potential, and p is the 
momentum operator. With 3C as the operator in (1), the eigenfunctions 
of 5C are the Bloch functions : 

Ko{v n (t; k) exp (tk-r)} = Sy(k)w*(r; k) exp (tk-r), (2) 

where 

p,(r; k) = u n (t; k) exp (tk-r), 

are the Bloch functions, and S*(k) are the energy eigenvalues. The 
subscript N refers to the different energy bands, and the function 
w w (r; k), which is normalized over the volume of the unit cell, has the 
periodicity of the lattice. In the presence of a static and dynamic field, 
the Hamiltonian becomes 
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= 2T„(P^ A ) !+F «- Fr > 



3c = 2 T^ p -; A ; + F(r) - F - r ' (3 > 

where e is the electronic charge, c is the vacuum light velocity, A is the 
vector potential for the dynamic field (the gauge <J> = and VA = 
is chosen), and F = cE is the force due to the static field E . The static 
and dynamic fields enter differently in (3), because of the fact that the 
vector potential for the static field increases linearly with time. It will 
be shown, however, that in the solution there is no distinction between 
the two fields. 

With the Hamiltonian given by (3), the wave function \p is chosen so 
that the equation of motion ( 1 ) does not contain secular terms, and so 
that the solution to the wave equation [(7) and (8) below] takes a 
relatively simple form in which the time-dependent and time-inde- 
pendent fields both appear. A normalized wave function that satisfies 
these two criteria is given by 

^ r >') = tt-^E [ d*lca N (t;k)u N [i;K(t)] 

(Z7T) 5 N J 

r i if i (4) 

X exptjkT + ^F-r* -- r J Q S»DE(r)]dr + Sjr(<;k)|, 

where A is the volume of a unit cell, and where, for the sake of sim- 
plicity, the function Sn(1; k) has been defined as 



S K (t',k) = it J drE-f dv u N *(r;K)V K u N (T;K). 



(5) 



Here V K is the gradient operator in K-space and the asterisk denotes the 
complex conjugate. One result of the added force on the electron, deter- 
mined by the total electric field vector 

!-*-!£, 

C dt 

is to change the electron wave vector from the constant value k to the 
time-varying value 

which describes the intraband electronic motion. The integral over k- 
space entering in (4) extends over all values of the wave number within 
an energy band, and the volume integral over coordinate space appearing 
in (5) extends over the volume of a unit cell. It should be noted that 
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the dc field does not appear symmetrically with the time-varying field 
in (4) ; this is due to the fact that the Hamiltonian has a different form 
for static and time-dependent fields. The form of the wave function 
given by (4) was suggested by Houston's choice of wave function for an 
applied dc field. 

It is assumed now that the wavelength of the applied radiation is 
much greater than the dimensions of the unit cell, and that the electron 
wave vector k has the same value for the lowest state in the conduction 
band as for the highest state in the valence band. Thus the electron 
momentum vector for the final state is the same as the momentum 
vector for the initial state (direct transitions). In Appendix A, it is 
shown that the substitution of the wave function \f/ from (4) into the 
wave equation ( 1 ) yields a differential equation for the undetermined 
probability coefficients a N (t; k). Considering only direct transitions 
between valence and conduction bands, there results [see (29)] 



dajdt = —Qa v 
da v /dt = Q*a c , 
where 



. e EM 

Q = i exp i 

m co cv 



I io cv dr — (S c — S v ) . 



(7) 



(8) 



The subscripts c and v refer, respectively, to the conduction and valence 
bands; M is the matrix element for transitions between the two bands, 



M = / dv Q u c *Vu v ; 



and io„ is an angular frequency given by 

S c - S u 

In taking the complex conjugate of Q it should be noted that the func- 
tion S N defined by (5) is real. This follows from the fact that the electric 
field E is chosen to be real, and that the integral 



/ 



dv u N *V K u N 
is imaginary, since 



V K / u N *u s dv a = V K 8 NN = 0. 
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The interaction between the semiconductor and the electromagnetic 
field results in the generation of a probability current density in the 
crystal. The expectation value of this current, in turn, acts as a source 
of radiation in the Maxwell equation for the curl of H. Therefore a 
complete solution of the problem of the interaction between the semi- 
conductor and radiation field requires an evaluation of this current 
source. The expectation value of the probability current density, J, 
can be expressed as follows, by noting that the canonical momentum 
operator p is Hermit ian, 



m 



7Y c*/ 



where V is the crystal volume and the angular brackets denote averaging 
over the volume V. 

It is now convenient to define a vector j by 



J = / d 3 k] 



such that j is the averaged current density per unit volume in k-space. 
In Appendix B it is shown that j satisfies the simple relationship [see 
(36)] 

j-E = dW/dt, (10) 

where W is the energy of the system per unit volume, i.e., 

W = (1/T)( |a c | 2 8 c + \a v \%). (11) 

Equation (10) is just the statement of conservation of power flow, for 
the left-hand side of this equation is the power transferred into the 
semiconductor, as given by the scalar product of current density and 
total electric field. The right-hand side of (10) is, of course, the resultant 
rate of increase of energy density in the quantized semiconductor sys- 
tem. 

Equations (7), (10), and (11) contain the necessary information to 
solve various types of semiconductor problems involving direct transi- 
tions. The derivation of these equations did not require discarding A 2 
terms in the Hamiltonian, and the solution for the unknown amplitudes 
a N did not necessitate using perturbation or approximation methods. It 
should be emphasized at this point that the results do not distinguish 
between static and dynamic fields. The set of equations (7), (10), and (11) 
is nonlinear and may be solved to any order in the electric field E to 
determine the various nonlinear interband phenomena. The remainder 
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of this paper will be devoted to a detailed treatment of the linear solu- 
tions, which will be shown to lead to the rate equation for the generation 
of electron-hole pairs (photoconductivity). 

III. SMALL-SIGNAL SOLUTIONS AND DERIVATION OF THE PHOTOCONDUC- 
TIVE RATE EQUATION 

The equations of motion for optical interband transitions [(7), (10), 
and (11)] can be greatly simplified by seeking solutions for a and a v 
correct to first order in the electric field E. Under this condition, which 
applies when the external electromagnetic force influencing the elec- 
tronic motion in the crystal is small, the wave vector K(f) can be re- 
placed by k, so that o) ev becomes time-independent. It is then easy to 
show from (7) and (8) that the solutions for a c and a v correct to first 
order are : 

a c = a c i0) -a w m -i?-M>[Ee i - l dt 

wrn. J (i2) 

o. = a v m - a c (0) -*L M*. ( Ee-^'dt. 

The superscript (0) indicates the zero-order solutions which, for this 
example, are independent of time. For the case where the crystal is in 
thermal equilibrium at 0°K, the constant amplitudes a o <0) and a„ (0> can 
be determined by noting that the electrons are confined to states in the 
valence band only, i.e., 

= 0, (13a) 



a, 



(0) 



and 

\a, m \Yk = p(k)cfk, (13b) 

where p(k) is the density of states. Integrating (13b) over a spherical 
shell gives for the number of states in a range die 

| a, w> | *d& = (V/ir 2 )k 2 dk, (13c) 

in which V is the crystal volume. 

A second-order linear differential equation can now be derived for 
the crystal current density j. Since the external forces are assumed to be 
small, K(0 » k, whereupon one obtains from (10) and (11) 

fro^ a) q c I 2 (14) 

1 V dt ' 
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where use has been made of the relation | a c | " -f- | a,. | " = constant. 
By substituting the first-order solution for a c from (12) into (14), it is 
found after some manipulation that j satisfies the equation of motion 
for a forced harmonic oscillator, i.e., 

S+-l-*£, as) 



where 



Ih e 



G= ■^V(-) \Me\ 2 \o. w \\ (16) 

and Me is the component of M in the direction of the electric held E. 

The validity of (15) may be checked by calculating the absorption 
coefficient a(w) for the fundamental absorption edge. It is assumed that 
the semiconductor has spherical energy surfaces with the same reduced 
effective mass, \/m* = \/m* — 1 /m*, along all crystal axes, and that 
in consequence 

We „ ^co„ + (ft/2m*)k\ (17) 

where fuaa is the minimum energy gap between the two bands. By making 
use of (13c) for | a„ <0> | 2 and (17) for w r „ , the observable crystal current 
density J can be found by integrating j, as determined from (15), over 
k-space. The resultant expression is given by the following constitutive 
relation between J and E 

where w is angular oscillation frequency of the radiation field. The 
absorption coefficient can now be found directly by observing that in 
the frequency domain the coefficient of E in (18) is the macroscopic 
conductivity, a, of the solid and is related to the absorption coefficient by 

a = Aira/llifi, 

in which ?< is the index of refraction of the crystal. The value for a 
obtained in this manner from ( 18) is the same as that given elsewhere. 7 
The photoconductive rate equation for generation of conduction 
electron-hole pairs also follows from the small-signal solution. If n is 
the number of conduction electrons per unit volume (equal to the 
number of holes), then one may express n as the k-space integral of 
l«c I 2 : 
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n = j d * k L?±l. (19) 

The rate of increase of conduction electrons is obtained from (14), 
whereby 



dn 
dt 



= E- (d'kj-. (20) 



In order to account for lattice collisions which produce intraband transi- 
tions and eventual interband recombinations, it is necessary to add a 
phenomenological loss term to (20). This is accomplished by adding a 
term n/r n to the left-hand side of (20), where t„ is the lifetime of an 
electron in the excited states of the conduction band. Equation (20) 
then leads to the photoconductive rate equation 



dn 



+ ^ = E- [d 3 k-i-. (21) 



A situation of particular interest is where E is an optical field made 
up of several coherent frequencies. Combining (15) and (21 ) then results 
in the rate equation for photomixing 

* + i.j«2:?M. (22) 

In this equation the subscripts p and q enumerate the different frequency 
components in the incident radiation. It is seen from (22) that n con- 
tains the sum and difference beat frequencies of the incoming radiation. 
This equation may be used to calculate the optical mixing properties of 
semiconductor crystals, where the mixing process results from photo- 
conductivity. 8 

IV. CONCLUSIONS 

The analysis of the problem of direct transitions in semiconductors 
has led to several simultaneous nonlinear differential equations. Lineari- 
zation of these equations correctly predicts the fundamental optical 
absorption and photoconductive mixing effects. Nonlinear effects such 
as might result from the time dependence of u cv ( Franz-Keldysh effect), 
or from the time dependence of | a c | 2 and | a v \ ''' (multiple-photon 
effects) can be predicted by appropriate higher-order approximation 
procedures. 
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APPENDIX A 



Derivation of the Equations of Motion Governing Direct Transitions 
Equation (2) can be rewritten as 

|- + 7(r)J u N (r;k) exp (tk-r) = S*(k)u„(r; k) exp (tk-r), (23) 

where Un(t; k) exp ikr are the stationary Bloch functions, and p is 
the operator — ifiV. It follows directly from (23) that in the presence of 
an electromagnetic field 

[24(p-c A ) 2+F(r) ]^ (r;K)eXP { k + ^0- r 

/ \ (24) 

= £ A .(K) W „(r;K) exp t/k + i FM ■ r. 

Equation (24) can be verified by performing the indicated operations 
and invoking (23). 

The substitution of the wave function ^(r, t) from (4) into the wave 
equation with the Hamiltonian of (3) yields, with the help of (24), 

- Z / A fe WArCr; K) + a„ \~- V K u N (r; K) 

" « (25) 

+ i ^ uAt; K) J| exp i(r, t, k, N), 

where exp i(r,t,k,N) is an abbreviation for the exponential function 

exp ijk-r + iF-ri - jj S. V (K) cfr + S N {t\ k) \ . 

Equation (25) is now multiplied by 

u Nf *(t; K') exp -i(r,t,k',N'), 

and integrated over all coordinate space. The resulting equation is then 
simplified by making use of the lemma 9 
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[ dv f d 3 kHr; k,k') exp i(k - k')-r = ^- [ cfc**(r; k,k), 

J 11 .p.o. J A J 

where $(r; k,k') is a lattice periodic function. Recall that the volume 
integral on the right-hand side of this expression extends over the 
volume of a unit cell. 
The result of performing the above operations on (25) is 



da 



} + ia L ^=-E«» #■{/ ******* 

', dt s ol {•> 



X exp i I ojln dr + (S N — Si,) , 



(20) 



(27) 



where Tuoln = Sl — Sat . It can be shown from the definition of S,. [see 
(5)] that the time derivative, dS L /dt, is 

— j = i— -• / dv UL*V K u L , 
dt at J 

so that (26) may be rewritten as 

^T= -Tt'av^-ifdvouSVKU*) 
dt N ot {>> ) 

X exp i J o) L k dr + (S N — S L ) , 

where J^' indicates that N = L is excluded from the sum. 
if 
There is a rather simple relationship between the matrix element of 
the gradient operator in K-space and in coordinate space, which for 
N 9* L, can be expressed as 9 

/ diw L *V K u N = — — / df u L *Vu N . (28) 

By combining (27) and (28) and noting from (6) that dK/dt = eE/fi, 
the results 



E-Y.' a * — { I dv Q u L *Vu N ) 

X exp i f u>,. N dr + (S N - S L ) \. 



da L __ . e 



For transitions between valence and conduction bands only, (29) reduces 
to (7). 
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APPENDIX B 



The Relationship between Current Density and Energy 

The expectation value, or volume average, of the probability current 
density in the presence of an electromagnetic field is given by 



j=^-:- a >- 



(30) 



Considering only transitions between the valence and conduction bands 
for the wave function of (4), the current density determined from (30) 
becomes 

J = — i — = / d 3 kliK(\ a c | 2 + I a v | 2 ) + | a c | 2 / dv u c *Vu c 
+ | a„ | 2 / dv Q u v *Vu v 4- a c *a„M exp i j co cv dr + (S v — S e ) (31) 

— a c a v *M* exp —i I a> cv dr + (S v — S e ) \>. 

Equation (31) can be simplified considerably by making use of the 
relationship 9 



i~ V*Stf(K) = iK 4- fdvou»*(r; K) Vu N (r;K). 



(32) 



Combining (31) and (32) results in the expression 

j = -i * f d'kii ™(|a. I 2 V A -S C + | a v | 2 V A S,) 
mV J { fi 2 

4- a c *a v M exp i f a cv dr + (S v - S c ) (33) 

— a c a„*M* exp — i / w c „ dr 4- (*S„ — S c ) \>. 

The energy density, W, for the quantized semiconductor system is 
given by (11). By differentiating W with respect to time, one has 



dl 



l y[\a c 



dK 

dt ' 



v*s c 



a v | — -V^S, 

at 



+ "^ + <"¥j 8 < + Kir 



* d fl « 



**>]■ 



(34) 
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Now, from the preceding work in Appendix A the time dependences of 
a c and a„ are known and, for the present situation, are given by (7). 
Thus by combining (7) with (34) and noting that dK/dt = eE/fi, 
one obtains the following formula for dW/dt: 

^ = ~|E-(| a c | 2 V*8 C + | a, | 2 V*S, 
dt Vn { 

— i — a c *a v M exp i I co c „ dr + (S v — S e ) (35) 

m Jo J 

h 2 T C l 11 

+ i— a c a v *M* exp — i I a> c „ dr + (S v — S c ) \f- 
m _h JJ 

A comparison of (33) and (36) illustrates the desired relationship 
between current density and energy density, viz., 

j-E = dW/dt, (36) 



where 



J - / d 3 kl 
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